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ABSTRACT 


We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves 
(LCs) display two humps at similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). 
SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, z = 0.0037, 20 Mpc). Optical and near-infrared photometry and spectroscopy are 
used to identify the energy source powering the LC. Nearly 50 epochs of high signal-to-noise-ratio spectroscopy were obtained within 130 days, 
o comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution 
<= of SN 2022xxf points to typical SN Ic/IcBL, with broad features (up to ~ 14000 km s7!) and a gradual transition from the photospheric to the 
+ nebular phase. However, narrow emission lines (corresponding to ~ 1000 — 2500 km s^!) are present from the time of the second rise, suggesting 
slower-moving circumstellar material (CSM). These lines are subtle, in comparison to the typical strong narrow lines of CSM-interacting SNe, for 
WJ example, Type In, Ibn, and Icn, but some are readily noticeable at late times such as in Mg I 45170 and [O I] 45577. Unusually, the near-infrared 
spectra show narrow line peaks, especially among features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. 
= We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario 
> is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and 
W) subsequent rapid LC drop. SN 2022xxf may be related to an emerging number of CSM-interacting SNe Ic, which show slow, peculiar LCs, blue 
CN colors, and subtle CSM interaction lines. The progenitor stars of these SNe likely experienced an episode of mass loss shortly prior to explosion 
© consisting of H/He-free material. 


——4 Key words. supernova: SN 2022xxf, ZTF22abnvurz 
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1. Introduction 2016] [Kangas et al.2017) Fang et al [2019} Pren- 
) that binaries play an important role in producing 
e e The demise of massive stars (Mzams = 8 — 10 Mo) as core- SESN S 
z collapse (CC) supernovae (SNe) comes in various flavors (see 
e.g. {Langer [2012] [Gal-Yam|2017). The diversity of SNe is In common SESNe, evidence for the presence of significant 
x. thought to be influenced mainly by the initial mass and the mass- circumstellar material (CSM) from progenitor mass loss is rare, 
C3 loss experienced by the progenitor star. Hydrogen-poor, stripped but has been found in some objects, e.g. late-time broad flat- 
envelope (SE) SNe originate from progenitors that have lost a topped Ha emission in a few Type IIb SNe 
significant part of their envelopes before the explosion. These in- Fremling et al.[2019). In the radio and 
clude SNe Type IIb (He-rich, little H), Ib (He-rich, no H), and Ic X-ray wavelengths, signatures of CSM are more frequently de- 
(no H and He). Significant mass loss, e.g. through strong stellar tected (e.g. . SNe 2014C (Milisavljevic et 
winds or interaction with a close binary companion, is required 2015), 2017ens (Chen et al.|2018), 2017dio (Kuncarayakti 
for a star to become a SESN progenitor. Such strong winds et al./2018), and 2018ijp (Tartaglia et al.[2021) constitute cases 
are expected from very massive progenitors (= 30 Mo, e.g. where Type Ib/c SESNe spectroscopically metamorphosed into 
[Groh et al.[2013), while in the binary scenario  CSM-interacting Type IIn supernovae, revealing the presence of 
the progenitors can be of relatively lower mass (x 20 Mo, e.g.,  H-rich external CSM. A small number of Type Ic SNe have been 
[Yoon|2015}[Dessart et al.|2020). Evidence is mounting from both i to nd with H/He-poor CSM, such as SNe oa 


studies of individual saan and samples (e.g./Taddia et al.|2015| 2014) and 20210cs (Kuncarayakti et al.]2022). 
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They show slow light curves with blue colors and distinct emis- 
sion lines due to the CSM interaction. (2020) 
presented two SESNe (20190ys and 2019tsf) that after a few 
months started to (re-)brighten. They concluded that the extra 
power needed for such a light curve (LC) evolution is presum- 
ably CSM interaction because none of the other powering mech- 
anisms at play at later phases are likely to result in such a be- 
havior. However, only SN 20190ys showed clear evidence for 
such interaction (e.g. narrow coronal lines) whereas SN 2019tsf 
did not. In the case of the bumpy LC of SN 2019tsf, interac- 
tion with a warped CSM disk influenced by a tertiary companion 


was suggested (Zenati et al.|2022). Interaction between a new- 


born neutron star and a binary companion star has also been pro- 
posed to produce bumpy SN LCs (Hirai & Podsiadlowski[2022). 
LC bumps are relatively common in (super)luminous Type Ic 
SNe. They have been suggested to be caused by central engine 
(e.g. magnetar) activities, CSM interaction, or combination of 


both (Gomez et al.|2021 | |Hosseinzadeh et al.|2022 
2022, Chen et al.|2023). 
In this paper we present SN 2022xxf, a SESN with a spec- 


tacular second hump! in the LC, but where the spectral se- 
quence does not exhibit strong CSM signatures. SN 2022xxf 
was discovered in NGC 3705 by [Itagaki] (2022) on 2022-10-17 
(MJDaiscovery = 59869.85) in a white-light image, and was re- 
ported to the Transient Name Server (TNS) on the same day. The 
host has a redshift of z = 0.00340 + 0.00001, distance modulus 
u = 31.54 + 0.45 mag, and luminosity distance 203. Mpc 
via NASA/IPAC Extragalactic Databasq?p. 
From measurements of host galaxy lines at the line of sight to- 
wards the SN (Sect. [2.2), we estimate and adopt a redshift of 
z = 0.0037 for the SN, which is used to correct the spectra. 
Spectral classification as a Type-IcBL SN was reported by [Bal] 
(2022), and confirmed with a spectrum taken earlier by our 


group (Nakaoka|2022} see Sect. 2.2). |Corsi et al.|(2022) reported 
a radio detection at 5.5 GHz using the Very Large Array (VLA). 


2. Observations and data reduction 


Observations of SN 2022xxf were done with a number of facil- 
ities, as listed in Table[A.1| which includes the instrument refer- 
ences. 


2.1. Photometry 


The first observations of SN 2022xxf/ZTF22abnvurz with the 


Zwicky Transient Facility (ZTF; Graham et al.|2019||Bellm et al. 
2019a), and the Palomar Schmidt 48-inch (P48) Samuel Oschin 


telescope under the twilight survey (Bellm et al.|2019b), were 
made on 2022-10-18 (MJDZTF = 59870.53) in r-band, 


first detection g . 
one day after the discovery. Photometry was obtained via the 
ZTF forced photometry servicd?] (Masci et al.[2019). No imme- 
diate pre-explosion non-detections are available as the SN just 


emerged from solar conjunction. After peak, additional g and i- 
band data were obtained with the ZTF camera on the P48 and 


! Tn the literature, the term “double-peaked LC" is used indiscrimi- 
nately for objects with early, fast-declining shock cooling emission e.g. 
SN 1993J (Richmond et al. [1994), or other types of slower peaks due to 
other mechanisms e.g. SN 2005bf RECETTE. Here we chose 
the word *hump' for SN 2022xxf as the rise and fall phases are well 
observed, forming roundish shapes, and the ease of association with the 
shape of a bactrian camel's back with the characteristic two humps. 

? 


ttps://ztfweb.ipac.caltech.edu/cgi-bin/ 


requestForcedPhotometry.cgi 
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Spectral Energy Distribution Machine (SEDM) Rainbow Cam- 
era on the Palomar 60-inch telescope. The P60 data were re- 
duced using FPipe for image subtraction. 
Further, we complemented the above with photometry from the 
Liverpool Telescope (LT), IO:O camera with ugriz filters. An 
automatic pipeline reduces the images, performing bias subtrac- 
tion, trimming of the overscan regions, and flat fielding. Tem- 
plate subtraction was done for the photometry. Photometry was 
also performed for the images taken with the 3.8-m Seimei tele- 
scope at the Okayama Observatory, Kyoto 
University, using the TriColor CMOS Camera and Spectrograph 
(TriCCS). Near-infrared (NIR) photometry was obtained using 
NOTCam at the 2.56-m Nordic Optical Telescope (NOT) at the 
Observatorio del Roque de los Muchachos on La Palma (Spain), 
and SOFI at the ESO New Technology Telescope (NTT) in La 
Silla, Chile. Standard reduction with bias subtraction and flat 
fielding was performed. No image subtraction was performed for 
the TriCCS and NIR data due to the lack of reference images. 


2.2. Spectroscopy 


The first spectrum of SN 2022xxf was obtained using the Hi- 
roshima One-shot Wide-field Polarimeter (HOWPol) on the 1.5- 
m Kanata telescope at the Higashi-Hiroshima Observatory, Hi- 
roshima University. Spectral classification as a SN IcBL using 
this spectrum was reported to the TNS (Nakaoka|2022). Within 
ZTF, a series of spectra were obtained with the SEDM and 
reduced with the pipeline described by (2019), 
and additionally with the Double Beam Spectrograph (DBSP) 
on the Palomar 200-in telescope, reduced using a DBSP re- 
duction pipeline (Mandigo-Stoba et al.12022) relying on Pypelt 
(Prochaska et al.|2020). Spectra were also obtained with the Al- 
hambra Faint Object Spectrograph and Camera (ALFOSC) using 
grism #4 on the NOT, by the ZTF and NUTS collaborations. 
ALFOSC data reductions were performed using ALFOSCGUP| 
We obtained optical spectra using EFOSC2, and near-infrared 
spectra using SOFI, both at ESO NTT, as part of the ePESSTO+ 
survey (Smartt et al.|2015). The raw data were reduced using 
the dedicated PESSTO data reduction pipelinqd?| Kyoto Okayama 
Optical Low-dispersion Spectrograph with optical-fiber Integral 
Field Unit (KOOLS-IFU) was also used for spectroscopy, using 
the VPH-blue grism. The data reduction was performed with the 
Hydra package in IRAF and dedicated softwarq'] All the spec- 
troscopic observations were accompanied with standard star ob- 
servations and followed by standard reductions which include 
bias and flat corrections, and wavelength and flux calibrations. 


In addition to low-resolution, we also obtained intermediate- 
resolution spectroscopy with X-Shooter at the ESO Very Large 
Telescope (VLT) in Cerro Paranal, Chile, on 2022-12-19] The 
observations were performed using the standard nod-on-slit 
mode, but each single arm spectrum was reduced using the 
“stare” mode reduction, given the brightness of the SN, and fi- 


nally stacked using the standard X-shooter pipeline (Goldoni et 
[al. [2006] Modigliani et al.|2010). Residual sky lines have been 


interpolated using the background as reference, and finally a tel- 
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luric correction was implemented using the line-by-line radiative 


transfer model (LBLRTM;|Clough et al.|1992). 


3. Results and discussion 
3.1. Extinction, metallicity, and redshift 


The amount of reddening in SN 2022xxf is estimated using the 
host galaxy narrow Na I D absorption lines in the intermediate- 
resolution X-Shooter spectrum. This yields Na I (D1+D2) equiv- 
alent width of ~ 1.5 A from the line measurements through 
Voigt profile fitting. Employing the relations from [Poznanski et] 
(2012), this equivalent width corresponds to an extinction of 
E(B— V) = 0.8+0.2 mag, assuming Ry = 3.1. We note, however, 
that in the range of the measured equivalent width the|Poznanski] 
[et al. |relation is not well sampled. Furthermore, any presence of 
CSM could have contributed to the Na I absorption and it could 
have a different Ry, complicating the extinction estimate. This 
caveat should be kept in mind where extinction is a factor in the 
analysis. The host galaxy of SN 2022xxf is inclined? at about 
68°, which makes high extinction likely in any case. The fore- 
ground Milky Way extinction is negligible, E(B — V) yw = 0.04 
mag (Schlafly & Finkbeiner|2011), which is confirmed by our 
measurements of the Na I D lines from the Milky Way in the X- 
Shooter spectrum. We thus assume E(B — V) = 0.8 + 0.2 mag as 
the total line of sight extinction for SN 2022xxf throughout the 
paper, unless mentioned otherwise. 

From detected interstellar medium lines of Ha and [N II] 
46584, the strong line method yields near-solar metallicity 
12+log(O/H) = 8.56 + 0.16 dex (N2, [Marino et al.|2013), al- 
though the strong presence of the SN spectrum prevents the de- 
termination of the host stellar population properties with cer- 
tainty. This metallicity is considered typical for a SN Ic, but high 
for a SN IcBL (Modjaz et al.[2020). Using these narrow lines, a 
redshift of z = 0.0037 + 0.0002 is derived and used to de-redshift 
the spectra. Note that this measurement corresponds to the ex- 
plosion site, which yields a slightly different value compared to 
the global host redshift value (z = 0.00340 + 0.00001; a differ- 
ence of ~ 90 km s7!). 


3.2. Light curves and color evolution 


The light curves of SN 2022xxf are presented in Fig. [I] Follow- 
ing the discovery, SN 2022xxf rises for more than 10 days to 
reach a maximum at r = 14.8 mag. While the very early rise 
is not well covered, the LCs around the first peak are consis- 
tent with those of typical SNe Ic/IcBL™| With g-band LC rela- 
tively fainter than the r and i bands, the color appears red, and 
the u-band detections are very faint. Assuming the host galaxy 
distance modulus (u = 31.54 mag) mentioned above and no ex- 
tinction, the SN magnitude at the first maximum corresponds to 
an absolute magnitude of M, — —16.8 mag. This is underlumi- 
nous for a SN IcBL, but still within the observed range for SNe 
Ic [Sollerman et l-2022). If reddening of 
E(B — V) = 0.8 mag is considered and corrected for (assuming 
Ry = 3.1), the LCs of SN 2022xxf become brighter, reaching 
absolute magnitude of ~ —19 mag, well within the range of SNe 


IcBL (e.g. Perley et al.12020). 


https://edd.ifa.hawaii.edu/inclinet/ 
10 Well-observed Typ gr (e.g. [Hunter et al.|[2009) 
Type IcBL SN 1998bw (e.g. 2001} |Clocchiatti et al. 


were chosen as representatives. The photometry of these and the other 
comparison objects were obtained from the Open Astronomy Catalog, 


https://github.com/astrocatalogs/OACAPI 


? Estimated using 


The color evolution of SN 2022xxf supports significant ex- 
tinction. As seen in Fig. [2| the dereddened (g — r) color of 
SN 2022xxf during the first hump is roughly consistent with that 
of typical SNe Ib/c around and after the main LC peak. If ex- 
tinction is not assumed, the intrinsic color of SN 2022xxf would 
be very red. An extra E(g — r) of ~ 0.3 mag would be required 
to bring the colors of SN 2022xxf to match the SN Ib/c tem- 
plate of (2015), although given the peculiar nature 
of SN 2022xxf it does not necessarily have to show the same 
colors as regular SNe Ib/c. 


In SN 2022xxf, following the first maximum (MJD 59880.0 
adopted as phase ¢ = 0 day), the LCs decline and at around +30 
days they rise again in all bands. The second rise lasts longer, for 
about 40 days, after which the SN reaches a second maximum 
at around +75 days. The second maximum is slightly brighter 
compared to the first one by 0.1 — 0.2 mag in the r and i bands, 
and clearly brighter (0.6 mag) in the g band. The u-band rise 
is dramatic (~0.1 mag d^), although there are only two data 
points, and with large error bars. This implies a color evolution 
towards the blue in the second LC hump, which is confirmed in 
the (g—r) color curve (Fig.[2). As the blue turnover is irrespective 
of the line-of-sight extinction, this is a robust observation. The 
peak brightness of the second peak is about 2 mag brighter than 
the luminosity expected from the ?9Co decay tail, equivalent to a 
more than six times increase of luminosity. 


To estimate the total radiative energy of the second hump, 
we calculate the bolometric LC of SN 2022xxf (Fig. using 
the method of [Lyman et al.) (2014). Extinction of E(B — V) = 
0.8 mag is assumed and corrected. Assuming that the first hump 
can be represented by a SN 1998bw-like LC, the bolometric LC 
of SN 2022xxf is subtracted by that of SN 1998bw, computed in 
the same way and scaled down to match the first maximum of 
SN 2022xxf (Fig.[A.1). The difference LC, which represents the 
second rise and hump, is then integrated along time, resulting 
in a total radiative energy of ~ 4 x 10° erg. This amounts to 
a significant fraction of the total radiative energy (^ 7 x 10% 
erg), but can be achieved in the interaction scenario that requires 
only a few percent of typical SN explosion energy of 10?! erg 
converted into radiation. Using the Hybrid Analytic Flux FittEr 
for Transients (HAFFET[!!] tool (Yang & Sollerman[2023). the 
first hump could be fit with a standard Ni-powered model with ~ 
0.4 Mo of CNi, while the slow rise of the second hump prevents 
a reliable fit with *°Ni. 


The second hump of the LC drops faster than it rises, which 
suggests that ?9Ni heating is unlikely to be the cause of the sec- 
ond hump; the decline is usually slower than the rise in the ?9Ni 
heating mechanism (e.g., see SNe 1998bw and 2007gr in Fig.[I). 
The production of ??Ni requires high temperature and density, 
thus this could only occur deep in the core, which implies that 
it is unlikely for the SN to synthesize any significant amount of 
?6Ni well after the explosion. Due to the longer photon diffusion 
time as the ejecta expand, >°Ni heating would predict a decline 
rate that is slower than the rise, which does not match the ob- 
servations of SN 2022xxf. Therefore, a different mechanism is 
likely to be at play. After the fall, the LC settles at the level of the 
tail luminosities of SNe 1998bw and 2007gr (after being scaled 
at the first peak luminosity), indicating that the first peak is pow- 
ered by the canonical *°Ni heating. There is however a hint of 
the LC flattening in SN 2022xxf (> +100 days), which suggests 
that extra radiative energy is still generated in SN 2022xxf in 
addition to the power from the radioactive decay input. 


!! https: //github.com/saberyoung/HAFFET 
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Compared with well-observed two-humped SESNe in the 
literature, SN 2005bf (Anupama et al. Tominaga et al. 
2006), its analogs PTFI Imnb (Taddia et 
2018) and SN 2019cad (Gutiérrez et al.]2021), the LC of 
SN 2022xxf appears to be distinct, although there is a diversity 
in the LCs of the other objects as well. In those, the second LC 
peak occurs 20 — 30 days after the first peak, i.e. in less than half 
the time compared to SN 2022xxf. Following the second max- 
imum, the LCs of the other objects decline slowly, whereas in 
the case of SN 2022xxf a more sudden drop in all bands fol- 
lows (~ 0.1 mag d^). As for the peak magnitudes, these other 
two-humped SNe peak at ~ —18 mag? at the second maximum, 
while SN 2022xxf is potentially brighter at ~ —19 mag, though 
we note the uncertainties in both distance estimates and in the 
extinction corrections. 

All these objects show a blue turnover in the color evolu- 
tion during the second rise, after which they become redder 
again. Only when the SNe approach the second hump do the 
colors become bluer than expected. This behavior may be inter- 
preted as a ‘normal’ radioactive first peak, followed by a sec- 
ond peak possibly powered by another mechanism generating 
extra energy. Several explanations for the powering of the sec- 
ond peak have been offered: an asymmetric explosion 


2006), a bimodal nickel distribution (Tominaga et al.|2005} 
Taddia et al.|2018 2021), or magnetar power 
(Maeda et al.|2007 . CSM interaction has 


not been invoked as a possible mechanism, due to the lack of 
strong narrow spectral lines, although it has been suggested that 
CSM interaction does not always require the production of nar- 


row emission features (e.g. 2001}|Sollerman et al.|2020| 


Dessart & Hillier|2022| |Maeda et al.|2023). The color evolu- 
tion of SN 2022xxf (Fig. |2) suggests that it becomes bluer as 


it approaches the second peak, starting from around +30 days. 
At the second peak, it has become bluer than the first peak by 
(g — r) ~ 0.7 mag, and subsequently the color becomes redder 
again as it falls from the peak. A similar blue turnover is seen in 
the other double-hump objects, although none of them reached 
colors bluer than (g — r) = 0 mag as in the case of SN 2022xxf, 
and they never became bluer than during the first peak. Such very 
blue colors are seen in SESNe interacting with H-poor CSM, 
such as SNe Ibn (e.g. 
2022). 
Ami et al.014j and SN 202 Tocs (Kuncarayakti et al. 022). 


3.3. Optical spectra 


Figs. [3] and [4] show the spectral evolution of SN 2022xxf dur- 
ing the first and second hump, respectively. With nearly 50 high 
signal-to-noise ratio spectra within 130 days, this dataset com- 
prises one of the most well-sampled observations of a (SE)SN to 
date. The spectrum initially appears smooth with broad profiles, 
typical for the early phase, and develops to show more features 
with time. The ejecta velocity at the pre-first maximum stage es- 
timated from the absorption minimum of O I 47774 is ~ 14000 
km s~!. Compared to some SNe IcBL such as SN 1998bw (e.g. 
Patat et al./2001), this is relatively low, but still well within the 
range of the SN IcBL sample (e.g. their 
Fig. 5). The strongest spectral features during the first hump 
are the P-Cygni profiles of Na I 445890,5896, Si II 46355, 
O I 47774, Ca II 448498,8542,8662 and the deep Fe absorp- 
tions at ~ 4200 and 5000 A, typically seen in Type Ic/ICBL SNe 


12 SN 2019cad could have reached nearly —20 mag if corrected for 


E(B - V) = 0.49 mag (Gutiérrez et al.[2021). 
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around peak brightness. This supports the notion that the first 
LC hump was powered similarly as the LCs seen in the majority 
of SESNe, via radioactive decay power and not by e.g. a shock 
cooling mechanism. During the second hump (Fig.H]. the spec- 
trum gradually morphs into becoming more nebular, as expected 


for SESNe a few months post-LC peak (e.g. [Patat et al.|2001] 
Hunter et al 009]. 


The overall spectral evolution of SN 2022xxf appears to be 
gradual, with no abrupt changes, and again similar to that of reg- 
ular SNe Ic/IcBL (see Fig. left). The transition from the 
photospheric to the nebular phase occurs relatively slowly, as in 
SNe IcBL, with the emergence of the nebular [O I] 446300,6364 
line at around +80 days. In the case of SNe Ic, the emergence 
of this line could occur earlier at around +60 days (Fig. 
or even before. There is no sign of H or He emission lines ap- 
pearing during the second hump, which is otherwise expected 
in the case of SN ejecta interacting with dense H/He-rich CSM 
(e.g. SNe oye SOT erar ORT aed SNe s EE Pe 
20190ys, and SNe Ibn, e.g. 
torello et al.|2007). When E(B — V) = 0.8 mag is considered, the 
dereddened spectra of SN 2022xxf appear to have an enhanced 
brightness at bluer wavelengths, relative to SNe Ic/IcBL at the 
corresponding epochs (Fig. left, +50 days onwards). This 
seems to support the CSM interaction interpretation. The rising 
blue continuum and strong Fe bump at around 5300 À are also 
frequently seen in the late-time spectra of interacting SNe, both 
in H-rich events such as SNe 2017dio and Type IIn SNe (e.g. 


SN 2020uem, {Uno et al.|[2023), and in H-poor events such as 
SNe Ibn (e.g. SN 2006jc, |Pastorello et al.|2007) and Icn (e.g. 
SN 2019hgp, |Gal- Yam et al.|2022). 


Looking at the spectra of the other two-humped objects 
(Fig. right), they and SN 2022xxf share some similarities 
during the respective LC phases, while clearly there are also dif- 
ferences. We note that the spectral classifications of these ob- 
jects are not identical: SN 2005bf was thought to be a Type 
Ic or Ib (or even possibly IIb with the interpretation of some 
features as hydrogen, [Anupama et al.|2005), PTF11mnb and 
SN 2019cad are both Type Ic SNe with lines narrower compared 
to those in SN 2022xxf. While the comparison is done for the 
same parts of the LC anatomy (first peak, valley, second peak, 
drop), the timescales are different for these objects. It is there- 
fore likely that the variations seen in the spectra reflect the dif- 
ferent time phases, and the emergence and disappearance of the 
second LC hump do not leave clear traces in the spectral evo- 


lution as previously pointed out (Folatelli et al.12006/ 
[al.[2018 2021). In those three objects, the LCs 


could be modeled with a bimodal Ni distribution (Orellana & 
2022) or with additional magnetar power input (Maeda 
2007 2021), although the effects of these 


mechanisms on the spectra are yet to be evaluated. 


3.3.1. Narrow emission lines 


In addition to the typical broad spectral features, a possible nar- 
row feature that is likely related to the second LC hump is iden- 
tified in SN 2022xxf. During the rise to the second peak (Fig. 4] 
from around day +46 to +73), the spectra of SN 2022xxf show 
a weak emission line at 7325 À, with a velocity FWHM of 
initially ~ 4500 km s^! which becomes narrower and reaches 
~ 2500 km s'!. This line then gradually disappears after the 
second peak, resulting in a flatter profile in that spectral region, 
before eventually the nebular [Ca II] 2A7292,7324 line starts to 
appear and take over (Fig. [A.3). The peak of the transient line at 
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Fig. 1: LCs of SN 2022xxf in ugrizJHK, bands. Apparent magnitudes (corrected for Milky Way extinction) are plotted in circles 
and squares, and dereddened absolute magnitudes (M, and M,, assuming E(B — V)tota = 0.8 mag and Ry = 3.1) in plus signs. The 
LCs of SNe 1998bw, 2007gr, 2005bf, 2019cad, and PTF11mnb are plotted for comparison, shifted in time to match their maximum 
epochs with the first maximum of SN 2022xxf, and in mag to match their first maxima. Vertical lines on top indicate spectroscopic 


epochs. 


7325 A is close to the 47324 component of the [Ca II] doublet, 
although this is not accompanied by the 47292 component. It is 
redder than the average wavelength of the [Ca II] doublet, and 
thus may be better associated with the [O II] 447320,7330 line. 
Narrow emission lines are more easily seen in the nebular 
phase of SN 2022xxf. Following the second maximum, narrow 
features emerge, superposed on the broad emission lines. All the 
broad nebular lines show a narrow core profile, indicating emis- 
sion from low velocities. This narrow line profile is seen in MgI 
45170, [O I] 45577, and all the major lines redwards i.e. Na I 
45893, [O I] 446300,6364, [Ca II] 447292,7324, O147774, OI 
48446, the Ca II triplet (all the individual triplet components), 
O I 49263, and other smaller structures across the spectrum (as 
shown in Fig. 5). These are unlikely to originate from the host 
galaxy. The narrow components are typically narrower than 2500 
km s^! (FWHM), superposed on broad components of ~ 5000 
km s~!. The peaks are generally offset by 300—500 km s^! from 
zero velocity inferred from the redshift, though this is compara- 
ble to the instrument resolutions in the case of the low-resolution 
spectroscopy. The narrowest emission lines such as Mg I 45170 
and [O I] 45577 display velocities up to +500 — 800 km s^! at 
their bases (HWZI, half-width at zero intensity), suggesting that 


they are unresolved given the instrument resolutiog!?] In addi- 
tion, other, mostly weaker, narrow lines are also found, including 
both known and unidentified lines, at 5532, 7010, 7155 ([Fe II]), 
7470, 7877/7896 (Mg II doublet), 8270, 8815, and 9436 (Mg D) 
À, some of which are seen already in the early nebular phase 
shortly after the LC drop, or even during the second rise in the 
case of 5532 À (Fig. [4). The [O III] 45007 line is variable and 
relatively broad, ~ +3000 km s^! at the base. In comparison, 
SNe Icn display narrow C/O emission lines corresponding to ve- 
locities 1000 - 2000 km s~! 
2022). and SNe Ibn show ~ 2000 km s™ in the He I emission 
lines (Pastorello et al.|2007). We do not detect high-ionization 
coronal lines (see e.g. 


Sollerman et al.[2020) in the intermediate-resolution X-Shooter 
spectrum, although it is likely that the CSM interaction was still 
weak at +50 days. The line profiles are generally Gaussian, with- 
out the wings of a Lorentzian profile, indicating that electron 
scattering is not significant. This is expected in ejecta dominated 


13 On the other hand, these velocities could be underestimated as the 
lines are situated on a pseudocontinuum, and thus their zero flux levels 
are uncertain. 
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Fig. 2: Color curve of SN 2022xxf in (g — r), considering the 
observed and dereddened cases, as compared to those of SNe 
2005bf, 2019cad, and PTF11mnb. Milky Way extinction is cor- 
rected for in all cases. The SN Ib/c color curve template from 


Taddia et al.|(2015) is also plotted for reference. 
p 


by intermediate mass elements, e.g. plasma of singly-ionized O 
will have 1 electron per 16 nucleon, therefore the electron scat- 
tering opacity in cm? g^! is 16 times lower than that of a H-rich 
plasma (1 electron per 1 proton). 


3.4. Near-infrared spectra 


The near-infrared (NIR) spectra of SN 2022xxf are presented in 
Fig. [6| Here, we follow the NIR line identifications of 
i. 2009); ho et al. (2021); Shahbandeh et aI. 2022). Relative 
to the optical spectrum, the NIR spectrum becomes nebular ear- 
lier due to the lower optical depth at longer wavelengths. In com- 
parison to other Type Ic/IcBL objects, while the rarity of NIR 
spectra of these objects especially at later phases only allows 
for limited comparisons, the earliest spectrum of SN 2022xxf at 
+51.8 days from X-Shooter displays a similar global appearance 
with broad emission lines near 10800, 11200, 11800, 13200, 
15000, and 16000 A. The main difference is the narrow peak 
of the Mg II 421369 emission line, which is not seen in the other 
objects, and a broad unidentified line at ~ 16850 À with similar 
strength as Si I 415888. The subsequent spectra show a similar 
set of features, as seen also in other SNe Ic/IcBL although pre- 
vious observations rarely reach later than +100 days. As in the 
optical, there are no clear detections of H and He lines. The fea- 
ture at 10800 À is attributed to a blend of C I/Mg II/O I and may 
contain He I 10830 Á, though not accompanied by He I 20581 
A, as in the case of SNe Ic/ICBL e.g. 
their Fig. 5). 

As the evolution progresses, the emission lines of 
SN 2022xxf gradually show narrower profiles (see Fig. [A.4). 
O I 411290 is clearly showing a sharp peak, and curiously a 
split line profile in the +79.7 days spectrum. Other lines also 
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show broad to narrow evolution, e.g. Mg I 411828 and 415033, 
as well as O I 413164. The broad profile at 10800 A initially 
shows a pronounced peak consistent with Mg II 210927, which 
slowly decays resulting in a flat-topped profile at later epochs. 
Following the LC drop (> +75 days), the Mg I lines become 
stronger relative to the O I lines. Mg II and C I also weaken dur- 
ing the same time period, which suggests cooling evidenced by 
the growing Mg I. Si I 415888 appears initially as a broad profile 
before developing a narrow peak, accompanied by an unidenti- 
fied stronger peak bluewards at 15770 A. The narrow features in 
the NIR spectra of SN 2022xxf are not seen in the other SNe at 
the corresponding epochs. They are coeval with those seen in the 
optical and show similar velocities, evolving from ~ 4000— 5000 
to ~2000 km s^!. 


4. Summary and conclusions 


SN 2022xxs displays an unprecedented LC evolution with two 
distinctive humps separated by ~ 75 days and a ~ 2 mag peak- 
to-valley amplitude, suggesting an extra energy input on the or- 
der of 4 x 10?? erg in addition to the regular *°Ni decay powering 
the first hump. The global optical/NIR spectral evolution is sim- 
ilar to the population of SNe Ic/IcBL, although the emergence 
of narrow features in SN 2022xxf during and after the second 
LC hump suggests the presence of a slower-moving material at 
1000 — 2000 km s~!. If this were due to CSM, the second hump 
may be explained by ejecta-CSM interaction producing extra ra- 
diative energy from the conversion of kinetic energy, although 
this did not affect the spectra significantly as the narrow emission 
lines are subtle. A CSM interaction scenario is also supported by 
the slow rise and fast drop of the LC in the second hump, and the 
flattening of the tail phase. 

While it is unlikely that the immediate vicinity of the pro- 
genitor star was completely free of CSM, and thus ejecta-CSM 
interaction could have taken place also in the early phases al- 
though relatively weakly, the second LC hump is naturally ex- 
plained by interaction with the bulk of the CSM. If the LC rise 
starting around 45 days after explosion (assuming that the explo- 
sion occurred ~ 15 days before the first maximum) corresponds 
to the inner edge of a detached CSM, then with a ~ 14000 km s! 
ejecta expansion velocity the location of the inner edge of the 
CSM would be around 5 x 10? cm. If such CSM were formed 
by material ejected from the progenitor star at ~ 2000 km s^!, 
this ejection must have occurred within ~ 1 year prior to the 
explosion. It is to be noted that the velocities seen in the nar- 
row lines may reflect a combination of expansion velocity of an 
unshocked CSM and the shock velocity in the ejecta-CSM inter- 
action, therefore both velocity components are in reality slower 
than the line width implies, which results in larger CSM distance 
and ejection time estimates. In the literature, it has been argued 
that mass-loss episodes within short pre-SN timescale could 
be caused e.g. by wave-driven outbursts (Wu & Fuller|2021), 
centrifugally-driven mass loss through spin-up (Aguilera-Dena 
et al./2018), or pair-instability pulsations (Renzo et al.[2020). 
Observationally, pre-SN outbursts on such timescales have been 
reported for the Type Ibn SNe 2006jc and 2019uo 
Strotjohann et al.[2021). For SN 2022xxf, we searched 


ZTF data up to 4.8 years prior to the SN with a total time cov- 
erage of 11% during this period, but no precursor eruption was 
found down to absolute magnitude of ~ —11 mag. 

SN 2022xxf may be related to the emerging subclass of SNe 
Ic with H/He-poor CSM interaction (‘Ic-CSM’), SNe 2010mb 


(Ben-Ami et al./2014) and 20210cs (Kuncarayakti et al.12022). 
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Fig. 3: Spectral sequence of SN 2022xxf during the first light-curve hump. Spectra are normalized by their average values and not 
corrected for reddening. Phases are in observer frame, relative to the first LC maximum. Prominent spectral lines are indicated with 
vertical lines corresponding to the rest wavelengths. Vertical grey shades indicate spectral regions affected by telluric absorption. 


These objects also show O and Mg lines in their spectra 
(Fig. [A.5). slow-evolving unusual LCs, and very blue colors. 
Their spectral evolutions are not well sampled, unfortunately, al- 
though SN 2010mb shows persistent nebular lines and a narrow 
[O I] 45577 emission line which is also seen in SN 2022xxf. The 
r-band LC of SN 2010mb shows a bumpy 180-days ‘plateau’ at 
M, ~ —18.3 mag, which could occur also in SN 2022xxf given 
the current 130-days evolution. It is, however, unclear how these 
SNe are related to the Type Icn SN class 
Pellegrino et al.|2022) although they all share the property of in- 
teraction with H/He-deficient CSM. The Ic-CSM objects show 
relatively regular photospheric to nebular spectral evolution and 
long-lasting LCs, in stark contrast to the current sample of SNe 
Icn. This suggests that the properties of the progenitors and the 
CSM could be very disparate. It is likely that the masses of 
ejecta and C/O/Mg CSM in SNe Icn and Ic-CSM are consid- 
erably different, where it could be low in the former and high in 
the latter (~ 3 Mo of CSM and z 10 Mo of ejecta in the case of 


SN 2010mb, [Ben-Ámi et al.[2014). 


The properties of the CSM and ejecta of SN 2022xxf are sub- 
ject to further study involving long-term monitoring and multi- 
wavelength observations (Izzo et al., in preparation). Observa- 
tions at later epochs may reveal additional clues on the origin 
and composition of the CSM, and therefore the associated mass 
loss of the progenitor star before the explosion. 
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Fig. 4: Spectral sequence of SN 2022xxf during the second light-curve hump. Spectra are normalized by their average values and 
not corrected for reddening. Phases are in observer frame, relative to the first LC maximum. Prominent emission lines are indicated 
with vertical lines corresponding to the rest wavelengths. Dashed vertical line indicates 7325 A, the average wavelength of the [O II] 
AAT320,7330 doublet (see sect.[3.3.1), and dotted horizontal line the approximate epoch of the second LC peak. Vertical grey shades 


indicate spectral regions affected by telluric absorption. 
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Fig. 5: Nebular spectrum of SN 2022xxf at +122.6 days. Insets show close-ups of the emission lines, with dashed vertical lines 


indicating the rest wavelengths of the emitting species. 
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Appendix A: Additional table and figures 
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Fig. A.1: Bolometric light curves of SN 2022xxf (red points) 
and SN 1998bw Clocchiatti et al.2011) (blue 
points; light blue for the scaled-down LC). Solid lines are spline 
representations of the LCs. The difference between the LC of 
SN 2022xxf and the scaled-down LC of SN 1998bw is plotted 
with dashed orange line. Dotted grey lines indicate ?Ni fits using 
HAFFET. 
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Table A.1: Facilities used in the observations. 


Telescope Instrument Band (A) R-A/AA Reference 

Kanata HOWPol 4500-9000 400 Li 
NOT ALFOSC 3500-9500 400 E R 

NTT EFOSC2 3600-9200 400 

NTT SOFI 9300-16400 600 

P60 SEDM 3800-9200 100 

P200 DBSP 3500-10500 1000 ORE X Gumi 1987] 
Seimei KOOLS-IFU 4100-8900 500 Maisubayashi et al. 
VLT X-Shooter 3000-24600 6000 Vernet et al.|(201T) 
Liverpool Tel. 10:0 ugriz — Steele et al. (2004) 

P48 ZTF Camera gri = E rar OTR 

P60 Rainbow Camera gri = 
Seimei TriCCS gri — ORL? 

NOT NOTCam JHKs — URL? 

NTT SOFI JHKs — (1998) 
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Fig. A.2: (Left panel) Spectral comparison of SN 2022xxf to SN Ic (green dashed lines) and GRB-less SN IcBL (blue dashed lines) 
templates (Modjaz et al.[2016), and other SNe after +50 days not covered by the templates. The template spectra are flattened, and 
thus do not represent the correct SED shape — comparison is intended for identifying similar spectral features. Well-observed Type 
Ic SNe 2007gr 2013ge are plotted in green solid lines, and SNe ICBL 1998bw 
and 1997ef (Modjaz et al.|2014) in blue solid lines. Spectra of SN 2022xxf are plotted in reddish colors, with the dereddened 
ones in lighter shades. Phases are in days after the light curve peak (first peak in the case of SN 2022xxf). (Right panels) Spectral 
comparison of SN 2022xxf with other two-humped objects, during specific LC phases: around the first peak, the ‘valley’ between 
the peaks, around the second peak, and the fall after the second peak. The LC-peak spectrum of SN 2007gr is shown in the first 
panel for comparison. 
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Fig. A.3: Evolution of the narrow line at 7325 À in velocity 
space. The rest wavelengths of [Ca II] and [O II] are indicated 
with dashed vertical lines (47324 is almost coincident with zero 
velocity). 
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Fig. A.4: Line profile of selected O and Mg emission lines in the optical and NIR. 
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Fig. A.5: Spectral comparison of SN 2022xxf to other interacting Type Ic SNe. 
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